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Abstract 

In the Randall-Sundrum model, a bulk neutrino field in the 5-dimensional space- 
time can give rise to tiny Dirac masses to neutrinos. In such a scenario, we have 
computed the contribution of the bulk neutrino field to the anomalous magnetic 
moment (g — 2)^ of muon. We have computed this contribution in the 't Hooft- 
Feynman gauge and have found that the contribution has the right sign to fit the 
current discrepancy between the experiment and the standard model value of (g — 
2)^. We have also studied possible constraints on the model parameters by including 
contributions to (g — 2)^ from other sources such as bulk gravitons. 



1 E-mail address: 
2 E-mail address: 
3 E-mail address: 



tprsh@iacs.res.in 

tpsr@iacs.res.in 

tpssg@iacs.res.in 



1 Introduction 



The models in extra dimensions have been proposed for solving the hierarchy between 
the electroweak and Planck scales [HE]- Among these the Randall- Sundrum (RS) model 
assumes the existence of two 3-branes connected by one extra spatial dimension and 
the metric in this model is non-factorizable [2J. The extra spatial coordinate has been 
orbifolded by the symmetry S 1 jZi and one of the 3-branes can be identified as the visible 
brane and the other as the Planck brane. The warp factor in this model suppresses any 
Planck scale quantities into electroweak scale on the visible brane. In the RS model, all 
the standard model fields are assumed to be confined on the visible brane and only gravity 
propagates in the bulk of the five dimensions. Subsequently, the RS model was generalized 
to include other bulk fields in order to explain physical quantities such as neutrino masses 
and mixing pattern which normally cannot be explained within the standard model of 
elementary particles. 

The deficit in the solar and atmospheric neutrino flux has given evidence for non-zero 
masses to neutrinos [3 J . Fitting to the data of solar and atmospheric neutrino experiments, 
the following set of mass-square differences are obtained: Am^ lai . = mj —m\ ph 7.6 x 10~ 5 
eV 2 and Am 2 tm = |m 2 — m\ 2 \ ~ 2.5 x 10 -3 eV 2 [I]. Here, m lj23 are the mass eigenvalues 
of the 3 active neutrinos. Apart from the mass-square differences an upper limit on the 
neutrino masses have been found through other experiments. Tritium /3-decay puts an 
upper limit on the neutrino mass scale to be of the order of 2 eV [5]. Whereas from the 
cosmological observations, the sum of the three neutrino masses needs to be less than 
about 1 eV [Sj. Since all the above experiments suggest tiny values for neutrino masses, 
perhaps a different mechanism should be operational for neutrino mass generation as 
compared to other fermion masses. 

To explain the smallness of neutrino masses in the framework of RS model, an addi- 
tional singlet neutrino field has been proposed, which like gravity field, propagates in the 
entire bulk of space-time [7]. The wave function of the bulk neutrino field is extended 
in the extra spatial dimension. The boundary conditions for this wave function can be 
chosen in such a way that it will have a very small overlap on our visible brane, resulting 
in tiny masses for neutrinos. The phenomenology of this model is determined through 
the Kaluza-Klein (KK) modes of the bulk neutrino field. Since the bulk neutrino field 
is singlet under the standard model gauge group, detection of KK modes of this field 
is challenging in the collider experiments. However, the loop effects due to these fields 
to any physically observable quantity can give us some hints about its existence. Here, 
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we study one of such observable quantities, namely the anomalous magnetic moment of 
muon. 

At the tree level, the anomalous magnetic moment of muon, ^-factor, has a value of 
2 and radiative contributions give some corrections to it. Hence, it can be quantified as 
a fi — (g — 2)^/2. For a review on (g — 2) M , see [91 [TO] . The world average value of (g — 2) M 
after the experiment E821 at the Brookhaven National Laboratory is as given below [5] 

ajf^ = 11659208.0(6.3) x HT 10 , (1) 

which is obtained with a precision of 0.54 parts per million. Various groups have computed 
the theoretical value for (g — 2)^ in the standard model. Most of the groups have found 
a discrepancy between the experiment and the corresponding standard model value of 
(g — 2) M at about 3a level [11]. Here, we take this difference as follows [10J 

Aa M = ajf p - a* M = (29 ± 9) x lO" 10 . (2) 

The above difference would indicate existence of new physics. Moreover, since this differ- 
ence is positive, the contribution due to new physics to (g — 2) M should yield a net positive 
value. 

As explained before, in the model of Ref. [7J, KK neutrinos can give some contribution 
to (g — 2) M , which we have computed in this work. Previously, some work in this direction 
has been done in Ref. [12] , where the authors have obtained a negative contribution to the 
(g — 2) M by adopting mass insertion approximation in the unitary gauge. In the present 
work, by carrying out an exact analysis in the mass eigenstate basis, we have computed 
the contribution from the bulk neutrino field to the (g — 2) M in the 't Hooft-Feynman 
gauge and scan the parameter space in the region consistent with experimental bounds. 
We have found that the bulk neutrino contribution to (g — 2) M has the right sign to fit the 
above mentioned discrepancy in (g — 2) M . We have also incorporated contributions from 
other sources such as gravitons to the (g — 2) M . Finally, we have studied the constraints 
that may arise from these various sources of (g — 2)^ in the parameter space defined 
in [7J. Recently, in [T3] . constraints have been obtained from lepton flavor violation by 
considering a set of models in the RS frame work. 

The paper is organized as follows. In the next section, we give a brief description of the 
warped extra dimensional model for neutrino masses [7J. In Sec. 3, we have computed 
the contribution of bulk neutrino field to (g — 2)^ in the 't Hooft-Feynman gauge. In 
this section, we have also performed a detailed phenomenological study of our obtained 
expression for (g—2) /1 . In Sec. 4, we describe possible constraints on the model parameters 
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of [7J by including contributions to (g — 2)^ from other sources such as gravitons. We 
conclude in Sec. 5. 



S= (fx dd)VG { E. 



2 Dirac neutrinos in the warped extra dimensional 
model 

As described before, the model [7J is based on the RS model, where an additional singlet 
neutrino field is introduced in the bulk of the 5-dimensional space-time. The metric in 
this model is 

ds 2 = e-^rj^dx^dx" + r 2 c d(f) 2 , (3) 

where <r(0) = kr c \(p\, r c is the compactification radius of the fifth dimension and k is an 
energy scale of the order of Planck scale, Mp. Here, <fi is an angular coordinate representing 
the fifth dimension. Due to the S 1 jZi orbifold symmetry, <p varies from to 7r and the 

3- branes located at these points are called Planck and visible branes, respectively. In this 
model, the invariant action for a singlet bulk fermion \l/ is as given below. 

- (^ a d A ^ - d A ^ a ^) + ^-^{ 7 a , a bc }^ - M fe sgn(0)M 
2 8 

(4) 

where G is the determinant of 5-dimensional metric Gab, is the inverse vierbein and 
bJbcA is the spin connection. The small case letters, a, b, c, run over flat 5-dimensions and 
upper case Roman letters run over curved 5-dimensional space-time. Here, M& is a bulk 
mass parameter which is O(Mp). The above action is invariant under — > —0, which 
should follow due to the Z 2 orbifolding of the RS model. This symmetry is known as 
0-parity which sets some boundary conditions on the wave functions of The bulk field 
^ decomposes into KK modes in the 4-dimensional world, which can be written as 

* L> kM) = J2^n R (x)^=ft R (<l>), (5) 
n v 

where ^l,r = |(1 T 7s)^5 V'n , ' R ( a ') are ^ ne ^ e ^" anc ^ right-handed KK modes in the 

4- dimensions and /„' (</>) are its corresponding wave functions. In order to have the 
following canonically normalized action for the KK modes in 4-dimensions 

«. = E/^{*.(»^*.w --».*.(»)*.(«)}. («) 

n J 

the wave functions f£ ,R should satisfy the following conditions 

dtfirwut) = f dtC(t)f«(t) = 6 mn , (±td t - v)ft R {t) = -x n tf*' L (t). (?) 



Here, m n are the masses of KK modes. In the above equations, the following change 
of variables have been done : t = ee a ^\ e = e~ kVc7T and — > y/kr c ef^' R (t). This 
change of variables imply that t = 1 corresponds to the visible brane. The unknown 
quantities in the above equation are : 

M b ™n krr7T /o x 

v = X' ^ = x e (8) 

The 0-parity, which is described above, imposes the following boundary conditions on the 
wave functions : fm*( e )fn( e ) = /m*(l)/?f (1) = 0- The coupled differential equations of 
Eq. ([7]) can be exactly solved and for n > 0, the wave functions can be expressed in the 
form of Bessel functions. The wave function for the zeroth mode can be written as 

fo L ' R (t) = fo L ' R W ±V , \fo' R (l)\ 2 = (9) 

The zeroth mode wave function of the right-handed field is highly suppressed on the 
visible brane, for v > ~. Because of this reason, by choosing the boundary condition as 
f^(e) = fniX) — 0, we can achieve a very small overlap of the singlet bulk field on the 
visible brane, for v > ~. Using this mechanism, in the next paragraph, we describe the 
neutrino masses in this model. 

The invariant action for neutrino Yukawa interaction is 

Sy = - J d 4 Xy/^~{%L (x)H {x)^ R (x, 7r) + h.C.}, (10) 

where L is a left-handed lepton doublet and H Q is the conjugate of the Higgs doublet. 
Here Y 5 has mass dimensions of -r= and the metric on the visible brane is (gvis)^ — 
g-M"")^ anc i g vis = M^vis)^. Both the lepton and Higgs doublets need to have the 
following rescaling : L = e^ a ^L, H = e a ^H, in order to have canonical kinetic terms 
for these fields. After substituting the KK mode expansion, eq. ([5]), the invariant action 
becomes 

Sr = -J d 4 xJ2vnLH^(x) +h.c, y n = VkY B f?(l) = Y B tf(l). (11) 

Here, Y 5 is an order one parameter. In the basis, ^ U L = {y L , ^f, • • • , \1>^) and = 
(^, #f , • • • , ^), we have 

( vy vy\ ■■■ vy n ^ 
mi ■•■ 

•-. 

m„ 



S Y = - I d 4 x^ u L M^ R + h.c, M 



(12) 



Here, the vacuum expectation value of the Higgs field is v — 174 GeV. Now, define 




(13) 



then the physical masses are given by 



L )^MV R = M t 



diag' 



(14) 



While diagonalizing the matrix M, the masses m n of KK modes are determined by the 
equation J v _i(x n ) = 0. However, the masses m n are not the physical masses of KK 
modes, since they get correction due to mixing with the left-handed lepton doublet after 
the electroweak symmetry breaking. After diagonalizing the matrix M, the lowest mass 
eigenvalue can be shown to be proportional to x/q which yields a very small value due to 
small overlap of the wave function f$ (1). The mass of light neutrino can be shown to be 



So far we have assumed the presence of only one bulk neutrino field and it generates 
one light Dirac neutrino mass. As explained in the previous section, to fit the neutrino 
oscillation data, we need at least two non-zero neutrino mass eigenstates. Hence, in the 
model of Ref. [7] we have to propose more than one bulk neutrino field. However, it 
has been argued in [7] that to cancel the anomalies related to 0-parity, only even number 
of singlet bulk fields can be introduced into the model. So, by introducing two bulk 
neutrino fields with slightly different bulk mass parameters, Mb, we can fit both the solar 
and atmospheric neutrino mass scales. In this picture the third neutrino has exactly 
zero mass. However, by introducing 4 bulk neutrino fields in an analogous way, we can 
generate all the three light neutrino masses. 

3 Contribution of KK neutrinos to (g — 2)^ 

We have shown in the previous section that to generate neutrino masses, the KK modes 
of bulk neutrino field have some mixing with the left-handed lepton doublet. This mixing 
would lead to gauge as well as Yukawa type interactions of muon with the physical KK 
modes of neutrinos. As a result of this, the KK modes of neutrinos give some contribu- 
tion to (g — 2) M , which we will describe shortly. It is to be noticed here that we have 
computed the contribution of bulk neutrino to (g — 2) M in the 't Hooft-Feynman gauge, 
hence the interactions of muon with Nambu-Goldstone states are necessary. For the sake 



[7\ 



m„ = 



V2u - lY 5 e {v -^v. 



(15) 
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of completeness, below we describe interaction terms involving the physical KK modes of 
bulk neutrino field. The interaction terms of other fields are unchanged from that of the 
standard model. 



The coupling of /j,W + uf k : 

In this model, the interaction term of VT-boson with charged lepton and neutrino would 
come from the 5-dimensional invariant action of the kinetic energy of the lepton doublet. 
For our particular case of muon field, these interactions are as given below. 



Sk3 j d*x^=[V^ 1 ^ L W+ +Wf^W- 
After using Eq. ( 1T5|) we get 



(16) 



S K 3 J d 4 x^= 



. .kt . u 1 75 



^.u 1 -75. 



(17) 



The coupling of fj,G + Uj k : 

Without loss of generality, we can go to a basis where charged lepton Yukawa couplings 
are diagonalized. In this basis, the necessary Lagrangian for our purpose is 



C = -Y^Hfi R -Y,ViWHVf + h.c. 



(18) 



Here L M is the left-handed lepton doublet containing muon. The second term in the above 
equation is from Eq. pip . We can use = m^/v where m M is the mass of muon. After 
using Eq. (IT5|) . we can show that 



c 3 -^Y/tGHut)"^+Y7FG + yt(K)"^ 



'7Ti7^+ (ttL\* 1 + 75 .. , Y^TTi7. 



* 1 - 75 



1,3 



a 2 



(19) 



3 h3 

The interaction terms, as described above, generate some contribution to (g — 2)^ at 
one loop level, which are shown in Fig. 1. The amplitude of any of the Fig. 1 is of the 
form iM. = —ieu(p')T tJ- u(p)e tl , where is the polarization of photon and e is the positron 
charge. After using the Gordon identity, u{p')^u{j>) = u(p') 



2m u 2m u 



q = p' — p, we can put the interesting part of the amplitude as 



1 W J 2m, 



-F 2 (q 2 ) + 



u(p), where 



(20) 



The contribution to (g — 2) M is Aa^ 



F 2 {0). 
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Figure 1: The contribution of KK neutrinos to (g — 2) /t . The momentum convention 
is that the outgoing muon has p', the incoming muon has p and the momentum on the 
photon line is q = p' — p, which is incoming. 



The contribution from Fig. 1(a) to (g — 2)^ is 



FJx) 



m, 



32tt< 



Xj), Xj 



w 



J=2 



1 



3 

31 19 11 2 7 3 fn . , 

1 x x H — x — (6 — x) mix) 

6 2 2 6 V ' K ' 



[m 



kk\2 ' 



(1-x) 4 

The sum of the contributions from Figs. 1(b) and 1(c) to (g — 2)^ is as follows 



F h (x) 



m, 



32n- 



X 



my? 



i=2 



' ~~ (mf) 2 ' 



(1 -xf 



— - + 2x — -x 2 — ln(x) 

_ _ 



(21) 



(22) 



S 



The contribution to (g — 2) M from Fig. 1(d) is 



AaJ?) = 



m. 



16tt 2 



E 



j=2,j=2,fc=2 



i 



kk 



F c (x) 
F d (x) 



1 

{l-xf 



[u^i (v*y + (u^y y k v£\ F d ( Xj ) 



11 2 1 3 i / \ 

- + —x — x + -x + x m{x) 

o A O 



m 



Xj — 



w 



(m^ k ) 2 ' 



1 1 2 , / \ 
x +xln(x) 

2 2 v ; 



(23) 



It is to be noticed that in the contribution to (g — 2) M from all the plots of Fig. 1, the 
summation in indices over k is from 2, since we have to subtract the light neutrino 
contribution which exists in the standard model. 

Here we comment on our results on the contribution of bulk neutrino to (g — 2) M . 
As already explained before that we have computed the contribution from Fig. 1 in 
the 't Hooft-Feynman gauge. The first three diagrams of Fig. 1 have given positive 
contributions, while the Fig. 1(d) has given both positive as well as negative contributions. 
In our numerical analysis, which we present below, we have found that in most of the 
parameter space the contribution from Fig. 1(d) is dominant and it gives a net positive 
contribution. Specifically, we have seen that the negative and positive contributions are 
comparable to each other but the magnitude of the former one is at least an 0(1) less 

(3) 

than the later one in AaJ, . Another comment is that, since there are two insertions of 
Yukawa couplings in Fig. 1(d) we would expect the contribution to (g — 2)^ to increase 
with Y 5 . In our numerical analysis we have found that this is true but there are some 
exceptions to this, which we will explain in our numerical results. 
The total contribution of one bulk neutrino field to the (g — 2) M is 



(24) 



The above contribution is mainly dependent on is, kr c and the 5-dimensional bulk mass 
parameter M^, which altogether determine the KK masses of the bulk neutrino field. 
Apart from this, Aa^ also depends on the dimensionless parameter Y 5 which determine 
the elements of the unitary matrices U L and V R . On the other hand, the light neutrino 
mass eigenvalue due to one bulk neutrino field, depends on is, kr c and Y§. As described 
in Sec. 1, from the neutrino oscillation data we have some idea on the magnitude of 
neutrino masses. For instance, in the hierarchical pattern of neutrinos at least two mass 
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eigenvalues should be : m so i ar = a/ Ar7ig 0lar w 0.01 eV and m atm = a/ Am atm w 0.05 eV. 
In this case, the third neutrino can have either zero mass or its mass should be less than 
of the order of m so i ar . In the case of degenerate neutrinos, all the three neutrinos should 
have a nearly equal mass and the common mass eigenvalue should be sufficiently larger 
than m atm . Since an upper bound from the cosmo logical observations indicate that the 
sum of the three neutrino masses should be less than about 1 eV [S] , we take the common 
mass to be as m^g ~ 0.3 eV. By fixing the mass eigenvalues of light neutrinos, we may 
eliminate v as independent variable. The parameter Y5 should be 0(1) and here we take 
its value in the range 0.1 to 2.5. Here the lower limit of 0.1 is due to the naturalness 
argument and the upper limit is due to perturbativity constraints. In order not to hit the 
Landau pole, the Yukawa couplings y n should be less than V^Lr. For n > 0, |/^(1)| = v2 
and hence Y5 should be less than y/2n 2.5. As for the bulk mass parameter M&, we take 
its value close to the Planck scale. We take the values of kr c to be around 12, since for this 
set of values we can produce TeV scale masses on the visible brane from the exponential 
warping. Finally, in our numerical analysis, we have allowed a 50x50 mixing mass matrix 
of Eq. fTT2l . After diagonalizing this matrix we get 49 KK neutrino mass modes which 
we have summed in the Aa^f . We have checked that the above mentioned number of 
KK modes in Aa^ is sufficient and its value vary insignificantly by further increasing the 
number of KK modes. 

In Fig. 2 we have plotted the contribution due to one bulk neutrino field to (g — 2) M . In 
Fig. 2(a) we have fixed kr c = 12.5. In this plot, the upper three curves are due to a bulk 
mass of Mb = 5 x 10 18 GeV and the remaining three lower curves are for Mb = 10 19 GeV. 
The three curves are arising depending on whether the light neutrino mass eigenvalue fits 
the solar, or atmospheric, or degenerate mass eigenvalue. The value of v is not shown in 
these plots, since as explained before it is not an independent parameter. The meaning of 
curves in Fig. 2(b) are same as that of Fig. 2(a), except that kr c has taken a value 12.0 
in Fig. 2(b). For one particular value of I5, the contribution due to a bulk neutrino field 
which fits the solar, or atmospheric, or degenerate neutrino is nearly the same. This is 
evident in both the Figs. 2(a) and 2(b), where the three curves due to different neutrino 
mass eigenvalues are close to each other. We can understand the reason for this as follows. 
For fixed values of Y5, kr c and Mb, the change in the value of Aa^ from one curve to 
the other curve can happen only due to the change in the neutrino mass eigenvalue. The 
ratio of neutrino mass eigenvalues between degenerate and solar cases could be at most 
by a factor of O(10). From Eq. ffT5l) . it can be seen that the neutrino mass eigenvalue 
is related to v dominantly through exponential factor, so we need only a change of ~0.1 
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Figure 2: Aa^ has been plotted as a variable of Y 5 . In these plots the three curves which 
are closely stacked together represents whether the light neutrino mass fits the solar, or 
atmospheric, or degenerate mass. The values of kr c and Mj, have been varied in both 
these plots. 

in v to produce a change of factor 0(10) in the mass eigenvalue. These exact numerical 
values in the change of v due to change in the neutrino mass eigenvalue can be seen in 
Tab. 1, where ranges of v are given for extreme values of Y5. Such a small change in v 
produces a slight change in the KK masses of neutrino fields, and hence a slight change 
to the Aa^. However, a closer examination would reveal that the amount of Aa^f where 
solar mass is fitted is slightly more than that for atmospheric mass which is even slightly 
more than that for degenerate mass. The change in these three cases could be at most 
in the first decimal place of Aa^ value. The values of v in these three cases would be 
different even though Y 5 could be same. We have presented the ranges of v for the cases 
of kr c = 12.5 and 12.0 in Tab. 1. The ranges of v in Tab. 1 are for extreme values of Y$ 





kr c = 12.5 


kr c = 12.0 


^atm 


v ~ (1.22 - 1.31) 
v ~ (1.18 - 1.26) 
v ~ (1.13 - 1.22) 


v ~ (1.25 - 1.34) 
v ~ (1.21 - 1.30) 
v ~ (1.16 - 1.25) 



Table 1: Ranges of v for different values of kr c and for different masses of light neutrino 
mass eigenvalues. In a particular range of the left- and right-end corresponds to Y5 = 0.1 
and Y 5 = 2.5, respectively. 
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of 0.1 and 2.5. Since the mass eigenvalue of neutrino, eq. f lT5|) . does not depend on the 
bulk mass parameter Mb, the range of v is also independent of this parameter. From the 
plots of Fig. 2, we have realized that there is some sensitivity between the light neutrino 
mass eigenvalue and the corresponding Aa^ value. However, for the same value of Y5, 
this sensitivity is so small that we may not determine the neutrino mass eigenvalue based 
on the Aa^f value. 

Next, we present results on how the Aa M varies by changing the kr c as well as the bulk 
mass parameter Mb, which is given in Fig. 3. In the plots of Fig. 3, we have fixed the 




Figure 3: Aa^ has been plotted against Y 5 for different values of kr c . In the top-left and 
top-right plots, the value of Mb has been taken as 10 19 GeV and 5 x 10 18 GeV, respectively. 
In the lower- middle plot, Mb = 10 18 GeV. The horizontal line in all these plots indicates 
the 2a lower limit of Aa^, see the text for details. 

light neutrino mass eigenvalue to be the atmospheric scale m atm . The choice of neutrino 
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mass scale do not make much difference in numerical values which we have argued around 
Fig. 2. From the plots of Fig. 3 we can understand that by increasing the value of 
kr c the exponential warping would decrease the KK masses of neutrinos and hence the 
contribution to Aa^ would increase. It can also be noticed from Fig. 3 that the lower 
the bulk mass parameter the larger the contribution is to the Aa^ , which is evident since 
the masses of KK neutrinos would become lower. For low values of Y 5 , the contribution 
to Aa 1 ^ is low and it is increasing and in some cases it may be saturated for large enough 
Y5. We have noticed that all the curves would be saturated for some large enough Y5. For 
example, the curves for kr c = 12.0 in Figs. 3(a) and 3(b) would saturate for Y 5 around 
50 and 20, respectively. The dependence of Aa 1 ^ on Y 5 is somewhat complicated as can 
be noticed from the theoretical expression given above. However, we can notice that Fig. 
1(d) has two insertions of Yukawa couplings and hence we can predict that Ad{, should 
increase with Y 5 . On the other hand, Y 5 also affect the contributions of the first three 
plots of Fig. 1 by determining the elements of U L and V R and also the physical KK 
masses of neutrinos. Numerically we have seen that for large value of Y 5 , the contribution 
from Fig. 1(d) is always 2 orders greater than that due to other plots of Fig. I. The 
contribution from Fig. 1(d) is given in the form of Aaj?\ from which we can see that 
there is a partial cancellation due to positive and negative contributions of Aa^K We 
have numerically seen that Aa^ goes to a saturation value for large enough Y 5 . Since 
in the plots of Fig. 3, we have fixed the neutrino mass eigenvalue to atmospheric scale, 
and the mass eigenvalue has an exponential dependence on u, we have found that after 
large enough Y" 5 the change in v would be far less compared to the change in Y 5 . Hence 
the values of m n which determine the KK masses of neutrinos would almost be saturated. 
The elements of U L and V R would change with F 5 , however, numerically we have seen 
that the net sum of the various KK modes is saturated after large enough value of Y 5 . In 
the plot of Fig. 3(c), for the case of kr c = 12.5 and = 10 18 GeV the amount of Aa 1 ^ 
would peak at around I5 ~ 0.1. In this particular case the contribution from Aa^ is 
significantly dominant at around Y5 ~ 0.1. However, in this case the lowest KK mass of 
neutrino is around 30 GeV. In other cases where the lowest KK mass is at least few 100 
GeV, Aa/j would give the dominant contribution. We have noticed these facts purely 
from numerical values. 

The results described above show that the contribution from a single bulk neutrino 
field in five dimensions can easily fit, depending on the values of parameters, the 2a 
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deviation of Aa M . From Eq. (jSJ) we take the 2a deviation as 

Aa^ = (1.1 -4.7) x 1(T 9 . (25) 

However, to be specific, the results of Fig. 3 indicate that the contribution to Aa M is most 
likely towards the lower end of the above 2a deviation or even less than this depending 
on the values of Y5. In the case of kr c = 11.5, the contribution to Aa^ is much below 
the 2a limit for M b = 10 19 GeV or 5 x 10 18 GeV. But keeping M b = 10 18 GeV and for 
Y5 ~ 2, we get a value of Aa M ~ 10 -9 . Since this is only from one bulk neutrino field and 
in a realistic scenario we need at least two bulk neutrino fields, so we can fit the above 
2a deviation by adjusting the M b and I5, even for a low value of kr c = 11.5. 

In the above analysis we have presented our results due to the existence of one bulk 
neutrino field. By introducing a second bulk neutrino field, the additional main parame- 
ters that the second field would carry are its 5-dimensional bulk mass parameter, M b , and 
its Yukawa coupling I5 to the muon. We can convince ourselves that the expression for 
Aa M due to this second field would be same as that of the first field, but replace the above 
said parameters accordingly. Now, consider a realistic scheme where there are two bulk 
neutrino fields which fits both the solar and atmospheric neutrino mass scales. In order to 
fit these neutrino mass scales, we may choose their Yukawa couplings to be nearly same 
but the values of v would be slightly different for these fields, which can be understood 
from Tab. 1. This would result in slight difference in the respective values of their bulk 
parameters M b . Hence the total contribution to (g — 2) M due to these two bulk fields 
would almost be the twice of the contribution from a single bulk field. 

4 Constraints on the model parameters by including 
graviton contribution 

In the previous section, we have shown that the contribution due to bulk neutrino field 
can fit the 2a deviation of Aa M . The fits in the previous section would in fact set limits on 
the parameter space of the model in [7], if this is the only source for (g — 2) M . However, the 
contribution due to the bulk neutrino may get further constraints due to the presence of 
other sources in the 5-dimensional warped model, such as from graviton or radion fields. 
Since gravity exists in the whole space-time, we cannot ignore its contribution to the Aa M 
[T4] . Similarly, the length of the fifth dimension should be dynamically generated and this 
leads to the presence of radion field. We have found that the contribution from gravitons 
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to Aa^ is significantly larger than that of the radion contribution [15]. Below, we show 
how the contribution from the gravitons to Aa M would lead to some constraints on the 
model parameters of [TJ. 

Regarding the graviton contribution to (g — 2)^, it has been computed in the case 
of RS model in [T6l [Hj. However, the expression given in [16J has been estimated by 
assuming all the standard model fields in the bulk of five dimensions. Since in the model 
of [TJ, only the singlet bulk neutrino and gravity are allowed in the bulk of space-time, 
we confine to the results given in [Hj, where it is claimed to be done for the case of the 
original RS model. In [14], an expression for the contribution of KK gravitons to the 
(g — 2)^ at one loop level is given, which is as follows 



where n c is the number of KK gravitons and A„. = e~ cW Mp. It has also been argued 
in [Hj that the number of KK gravitons are bounded by some unitarity constraints due 
to 77 elastic scattering, where it has been shown to be n c < 10 — 100. By taking the 
Planck scale as Mp = 10 19 GeV, the contribution due to a single KK graviton are as 
follows : (kr c , Aa£) = (12.2, 6.8 x 10~ 9 ), (12.1, 3.64x 10" 9 ), (12.0, 1.94 x 10" 9 ), (11.7, 2.95 x 
10~ 10 ), (11.3,2.39 x 10 -11 ). Comparing these values with the 2a range of Eq. (j25j) . 
anything above the kr c = 12.1 can be ruled out purely from the KK graviton contribution 
to {g— 2)^. For kr c = 12.0, only two KK gravitons can exists. Whereas, for kr c = 11.7 and 
11.3, of the order of 10 and 100 number of KK gravitons can exist, respectively. Although 
it may be inappropriate to assume the existence of only one or two KK gravitons in the 
universe, nevertheless, by the unitarity bounds of [14J it seems to be consistent. In these 
cases we show how this assumption would put bounds on the bulk neutrino parameters 



For kr c = 12.1, there can exist only one KK graviton to fit the 2a deviation of Aa M . In 
this case, an amount of Aa™ ax — Aa^ = (4.7 - 3.64) xl0~ 9 = 1.06xl0~ 9 can be shared by 
contributions due to other sources. It has been shown in [15] that for radion mass greater 
than about 200 GeV, its contribution to Aa M would be less than 10~ 9 . For simplicity, 
we ignore this contribution by assuming radion mass to be greater than 200 GeV. Then 
in the degenerate masses of light neutrinos, since at least 4 bulk neutrinos should exist, 
the contribution from any single bulk neutrino should be Aa^ < 2.65 x 10~ 10 . On the 
other hand, if the neutrino mass pattern is hierarchical then at least two bulk neutrinos 
should exist. By assuming that each bulk neutrino contributes by the same amount to 
(g — 2) M , we can then put a bound Aa^ < 5.3 x 10~ 10 . An upper bound on Aa^ would 




(26) 



of 0. 
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Deg 


enerate neutrinos 


Hierarchical neutrinos 


M b 


Y 5 


V 


Va 
I 5 




ys 


V s 


10 19 GeV 


1.2 


1.22 


1.9 


1.28 


1.8 


1.32 


5 x 10 18 GeV 


0.6 


1.2 


0.9 


1.26 


0.9 


1.31 


10 18 GeV 


0.1 


1.15 


0.14 


1.21 


0.14 


1.25 



Table 2: Upper limits on Y" 5 and v for different values of M& and in different patterns of 
neutrino masses. In the case of hierarchical neutrinos, the suffix a and s represent for 
atmospheric and solar neutrinos, respectively. 

imply an upper limit on Y 5 and thereby an upper bound on the corresponding v of the 
bulk neutrino. In Tab. 2 we show conservative bounds on these parameters for various 
values of bulk parameters In the case of kr c = 12.0, two KK gravitons can exist to 
fit the (g — 2)^. The above analysis done for kr c = 12.1 can be repeated for kr c = 12.0 
and we can get bounds on the parameters Y 5 and v accordingly. However, if we include 
the contribution of the order of 10 KK gravitons then kr c should be close to 11.7. In this 
case, by comparing the curves for kr c = 11.5 of Fig. 3, for Mb close to 10 18 GeV, a single 
bulk neutrino can contribute an order of 10~ 9 to (g — 2) M . So for this choice of Mj, we 
can get upper limits on Y 5 and v. However, if the number of KK gravitons is of the order 
of 100 then kr c should be around 11.3. For this low value of kr c , a single bulk neutrino 
contribution would be at most 6 x 10 -10 even for Mb = 10 18 GeV. Hence the constraints 
on the bulk neutrino parameters would be less stringent in this case. Finally, to comment 
on the radion contribution to (g — 2) M , the constraints on the bulk neutrinos which are 
described above may become even severe if the radion mass is less than about 200 GeV. 

In the previous paragraph we have described constraints on the bulk neutrino param- 
eters which are arising purely from the 2a deviation of Aa M . However, these constraints 
can be even stringent by including other observable quantities such as BR(fi — > ej). 
The current experiments have not found the decay fi — > ey and put an upper limit on 
its branching ratio as BR(fi — > ej) < 2.4 x 10~ 12 at 90% CL [17] . In the context of 
the warped model [7], this decay channel has been studied in [TS]. To satisfy the upper 
bound on BR(fi — > ej) and to get an appreciable contribution from Aa^, we may have 
to fine-tune the Yukawa couplings. It is interesting to study these effects but this is out 
of the scope of the current work. Finally, we comment that in this work we have studied 
correlation between neutrino mass eigenstate and (g — 2) M for a single generation. We 
hope to extend this to include all the generations in a future work. Some related work in 
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the context of different models has been studied in [T9] . 

5 Conclusions 

The RS model, which is based on the warp geometry, is elegant to explain the hierarchy 
between the Planck and electroweak scales [2]. Subsequently it was shown that such 
model can offer possible explanation of the origin of small neutrino masses when the right 
handed component of the neutrino is allowed to propagate in the bulk [7] . The signals of 
such a model lies in the detection of KK modes of the bulk neutrino field. In this work, 
we have studied an indirect signal of these KK modes, i.e. the virtual effects of these KK 
modes to the anomalous magnetic moment of muon, (g — 2) M . We have computed the 
contribution of bulk neutrino to (g — 2)^ at one loop level in the 't Hooft-Feynman gauge. 
Comparing our result with the experimental limits, we have found that this contribution 
comes with the right sign so as to fit with experimental data. Moreover, in some of the 
parameter space of the model in [7], the 2a discrepancy of (g — 2)^ can be accommodated. 
We have also studied how the bulk neutrino contribution would be restricted by including 
the graviton contribution to (g — 2) M . 
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